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Abstract 

We perform a study of the final state interactions of the K^K^ and the K^d systems in the 
reactions pp — > dK^K^ and pp — > dTr~^r]. Since the two-meson system couples strongly to 
the ao(980) resonance, these reactions are expected to be an additional source of information 
about the controversial scalar sector. We also show that these reactions present peculiar 
features which can shed additional light on the much debated meson-baryon scalar sector 
with strangeness —1. We deduce the general structure of the amplitudes close to the dK~^K^ 
threshold, allowing for primary K^K^ as well as tt+t/ production with the two mesons in 
relative S- or P-wave. The interactions of the mesons are accounted for by using chiral uni- 
tary techniques, which generate dynamically the ao(980) resonance, and the K^d interaction 
is also taken into account. General formulae are derived that allow to incorporate the final 
state interactions in these systems for any model of the production mechanism. We illus- 
trate this approach by considering two specific production mechanisms based on three flavor 
meson-baryon chiral perturbation theory. It is demonstrated that in this scenario the K^d 
interactions are very important and can change the cross section by as much as one order of 
magnitude. The amount of 7r"'"r/ versus K^K^ production is shown to depend critically on 
the primary mixture of the two mechanisms, with large interference effects due to final state 
interactions. These effects are also shown to occur in the event distributions of invariant 
masses which are drastically modified by the final state interactions of the two-meson or the 
Kd system. 
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1 Introduction 



The reaction pp — > dK~^K^ is presently the subject of experimental study by the ANKE collabo- 
ration at the Cooler Synchrotron COSY at Jiilich with the aim (among others) of learning about 
the nature and properties of the ao(980) resonance [Q. The problem has attracted also the interest 
of theoretical groups @,§(see furthermore the contributed papers in [Q). The prospect of gaining 
novel information about the ao(980) resonance, which might help to shed further light from the 
experimental side on the disputed nature of this resonance, is one of the attractive features of this 
reaction. This controversy originates from the observation that there are several different models 
to deal with the isospin / = 0, 1 scalar sector, all of them reproducing the scattering data to 
some extent, but with different conclusions with respect to the origin of the underlying dynamics. 
In particular, in refs. 1^-0 these resonances are considered as preexisting ones (genuine quark 
model states), while in ref. |jlO| they appear as meson-meson resonances generated by a potential. 
In ref. the advocated solution is that the ao(980) and the /o(980), the latter in the channel 
with 1 = 0, are exotic resonances, that is, not simply qq states, while the preexisting qq scalar 
nonet should be heavier, around 1.4 GeV or so. Other interesting approaches to this problem 
are refs. [lT2|-p!^, the relativistic quark model with an instanton-induced interaction of the Bonn 
group |T^, the Jiilich meson-exchange approach |]TB[, the Inverse Amplitude Method |T^ or some 



variants of it EO]. It is notorious that opposite conclusions are obtained in refs. [pT|,[2^ from the 



use of QCD sum rules. Regarding this controversy about the nature of the scalar resonances, the 
works of refs. [23-28] have stressed the role of chiral symmetry and unitarity to understand the 
dynamics associated with the lowest lying scalar resonances (see also ll^). As a result of the 
latter references the lightest O"*"^ nonet is established to be of dynamical origin, i.e. made up of 
meson-meson resonances, and is formed by the (t(500), k, ao(980) and a strong contribution to the 
physical /o(980). On the other hand, the preexisting scalar nonet would be made up by an octet 
around 1.4 GeV and a singlet contributing to the physical /o(980) resonance. Similar conclusions 
about the nature of the scalar resonances have been obtained in ref. PT| for the channels with 
isospin and 1 within QCD sum rules. The previous set of works [p3| -p8| constitute a unique 
theoretical approach to the scalar sector capable to study all these reactions in an unified way. 
This is achieved because all these processes are related by the use of an effective theory description 
that combines chiral perturbation theory and unitarity constraints. 
It is important to remark, as already pointed out in ref. 



12], that it is mandatory to study 



not only the experimental data concerning phase shifts and inelasticities related to meson-meson 
scattering but also production reactions where the interactions between the mesons manifest 
themselves through the final state interactions (FSI). This was indeed the subject of refs. lp5|-[28[] 
and plays a key role in the present investigation. 

In the baryon sector the studies of )p9|-|31[] based on the use of chiral Lagrangians also show that 
the A(1405) resonance is generated in a similar way. The fact that the A(1405) is of dynamical 
origin was pointed out already many years ago ]]32|. More recently it has also been shown iST 



that the A(1670) and the S(1620) are generated dynamically in the same chiral scheme of ]|30 ]. 
Still, in this sector more work is needed to firmly establish these results. 

As we will see, the reaction pp dPQ (where P, Q denote pseudoscalar mesons) offers novel 
possibilities with respect to other reactions where the ao(980) resonance is produced since it is 
sensitive to both the meson-meson and meson-baryon final state interactions. In particular, we 
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stress the importance of the Kd FSI, which has a very pronounced influence on observable yields, 
invariant mass distributions or cross sections, mostly through the interference with the meson- 
meson FSI (which gives rise to the ao(980)). This should provide extra information to test the 
implications of chiral symmetry on the nature of the low-lying resonances. 

The manuscript is organized as follows. In section H, we discuss the basic reaction mechanisms 
for the process pp dPQ, first in very general terms and then we consider a specific model based 
on chiral symmetry for the primary production of the meson pair. The final state interactions are 
treated in section |^, separately for the meson-meson and the meson-baryon systems. We stress 
in particular the role of the anti-kaon-deuteron FSI. The results are presented and discussed in 
section ^ and conclusions are drawn in section ^. The appendix contains a detailed discussion on 
the general structure of the process pp dPQ. 

2 Basic reaction mechanisms 

2.1 General considerations 

The reaction measured in |1| is: 

pp dK^K^ . (1) 
We will study it theoretically in connection with the accompanying process: 

pp — >• dr^^rj , (2) 

since the dynamics of coupled channels, which we will use here, deals with both channels simulta- 
neously. On the other hand, the energy of the ANKE experiment is fixed to ^/s = 2912.88 MeV 
just about 45 MeV above the dK~^K^ threshold. 

The reaction (P forces the K^K^ system to be in an / = 1 state which, given the proximity 
of the ao(980) resonance, would have its rate of production and invariant mass distributions 
very much influenced by the tail of that resonance. The reaction (^, which is also planned to 
be measured by the ANKE collaboration, could see the actual shape of the ao(980) resonance 
through the mass distribution of the vr"*"?] system. 

The P-wave nature of the reaction [Q, 0] is another peculiar feature that makes it different to 
other ones producing the ao(980) |^,|3^. Indeed, due to total angular momentum and parity 



conservation as well as to the antisymmetry of the initial state, the two mesons cannot be simul- 
taneously in intrinsic S-wave and in S-wave relative to the deuteron. Henceforth, we denote by i 
the orbital angular momentum of the CM motion of the two pseudoscalars PQ and the deuteron, 
and by L the orbital angular momentum of the pseudoscalar mesons in their own CM frame, what 
we also call intrinsic angular momentum of the two mesons. With this notation the cases i = 1, 
L = and i = 0, L = 1 are possible and then the initial state is forced to have io = 1, 3 and 
S = 1, with io denoting the orbital angular momentum of the initial pp state. These are the 
dominant contributions since they imply the lowest power, namely 1, of the small three- momenta 
of the deuteron or kaons in the transition amplitudes. Of course, the threshold of the reaction 
is much lower than the one of ([^) but due to the resonant nature of the interactions between 
the pseudoscalars, which only occurs when both 7r+?7 and K'^K^ are coupled together |[T9|,p3|,|2^, 
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we will consider the same structures as well for (|^). Even more, since the tt+t] system is expected 
to couple only very weakly to the P-waves then only the structure with £ = 1, L = is kept for 
TT'^ri. The fact that either £ = 1 or L = 1 has its relevance since, close to threshold, one would 
expect to have all particles in an S-wave. Thus, the K~^K^ system would be in an intrinsic S-wave 
state subject to the full strength of the ao(980) resonance tail. However, in the present case the 
contributions with L = or 1 are of the same order of magnitude. Furthermore, as commented 
above, if the tt^i] system does couple just weakly to the P-waves, as we will argue below that 
this is the expected behaviour, then the reaction mechanism establishes a distinction between the 
K~^K^ and Tr~^ri production processes which is novel with respect to other threshold production 
reactions. 

The reaction mechanisms for K^K^ and tt+t] production are complicated. Some approximate 
models have been used in refs. based on direct production of the ao(980) resonance. The 
7r+7y production channels are not studied there although they could in principle be accounted 



for by using partial decay rates of the ao(980) into KK and vrr^ |^6|]. However, the arguments 
given above about the unique prospects of this reaction indicate that the rates could be quite 
different in the present process than those determined by making use of standard Breit-Wigner 
parameterizations or related ones. Thus, substantial deviations from the rates observed in the 
present reaction would give further support to the method of coupled channel generation of that 
resonance as the appropriate tool to deal with the light scalar resonances, versus static pictures 
that consider this resonance as an object with predetermined decay rates into different channels, 
and the final state interactions are just taken into account via Breit-Wigner like modifications of 
the resonance propagator. Indeed, inconsistencies in the treatment of the related /o(980) resonance 
as a pure Breit-Wigner or following a Flatte formula have been recently pointed out in ref. ^\ in 
the reaction 'yir^iT^. 

2.2 Dynamical model based on chiral symmetry 

Since we are interested in stressing the role of the final state interactions of the KK and Kd systems 
we refrain from searching for a complete model and simply parameterize the original structure 
of the amplitude close to the dKK threshold, the region of interest for the ANKE collaboration. 
This is done in Appendix A. Nevertheless the number of free parameters when allowing for the 
general structure deduced in Appendix A is too large to draw any definite conclusion. To overcome 
this difficulty we consider in the following a specific model derived from lowest order three fiavor 
chiral perturbation theory (CHPT) |]38|- ^] depicted in fig.|I|, giving rise to two definite structures 



which already illustrate the main points of our investigation, namely, the extreme importance of 
the final state interactions between the mesons and between the K^ and the deuteron. Then a 
variety of observables are evaluated in terms of the two parameters of the theory, up to a global 
normalization which will not be needed to evaluate ratios of cross sections and invariant mass 
distributions nor to compare with the future ANKE results. It is worth stressing that, given any 
model accounting for the primary production of the dK^K^ and (ivr+r/ systems, the final state 
interactions can be taken into account following the general scheme presented in Appendix A. 
Let us come back to fig.|l|. Such diagrams involving the production of two mesons were evaluated 



in ref. |^T| and have a similar structure, with some cancellations among them. In order to see this 



structure we consider diagrams |I]c) and|I|e). The two-baryon-three-meson (BBMMM) vertices are 
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given by ^ 



C[^^ = \{D + F) {nYl.u^^'^P + Pl'l.u^i'h) , (3) 



where is a SU(3) matrix containing the meson fields which is given exphcitly in [^. Fur- 
thermore, D ~ 3/4 and F ~ 1/2 are the canonical SU(3) axial coupling constants. The relevant 
terms for our case are contained in the w^f and u]^^^ matrix elements, and keeping in mind the 
non-relativistic reduction of 7^75, are proportional to: 

wf) - ^7^'{d,K'K- -d,K-K'^)+VQd,r^K'K- , 
v2 

d^Tv+K'^K- + Ti+d^K-K"^ - 2TT+d^K'^K- . (4) 
Thus, diagram |l]c), with a non-relativistic reduction of 7'^75 leads to a structure 

-^o^^\pk+ -PRa) , (5) 

which produces the KK system with a relative P-wave, while diagram |l|e) leads to a structure 

V6a(2)g, (6) 

where the superscript 1(2) applies to the the baryon line to the left (right) side of fig.|l]. In addition, 
the initial proton on the left (right) baryon line has three-momenta Pi{p2)- 

The global structure of the amplitude is obtained by considering also the a'^^^'q vertex in the 
meson baryon vertex to the left side of the diagram. By taking q = pi — pd/2, with pi, pd the 
momenta of the initial proton and the deuteron, respectively (we have checked that consideration 
of the Fermi motion in the deuteron does not change the final structure), we find two types of 
terms: 

I) a(i)(pi - pj2) a^^\pK+ - pro) - ff^'\pK+ - Pro) ff^'\p2 - Pd/2) , 

II) (pi - pd/2) (pi - pj2) - {p, - Pd^ ^^'^ {P2 - pjl) , (7) 

which would come from the sum of the diagrams figs.|I]c) and |I]f), for eq.(^I), and diagrams 
figs]l|e) and Pi), for eq.(^II), after taking into account the isospin zero of the deuteron. Also the 
deuteron wave function appears with its value at the origin, 0(0), in coordinate space neglecting 
the range of the interaction, since the propagators of the pseudoscalar mesons in fig.|^ carry very 
high momentum transfers. In order to take into account the antisymmetry of the initial pp state 
we must subtract to the previous expressions the same amplitudes exchanging the two initial 
protons. This means both spin and momentum but, since we have = 1 in the initial state, 
the wave function is spin symmetric and then it is enough to subtract the amplitudes of eq.(|^ 
exchanging pi P2 = — Pi- This leads, up to a global factor two, to the structures: 

I) a^^\pK+ -Pko)+ ^^^\pk+ - pro) ^^'^Pi , 

II) . (8) 
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This calculation would not be complete to account for the vr exchange because the explicit 
use of the isospin deuteron wave function forces also the simultaneous consideration of diagrams 
figs]l|d) and |l]g) with the exchange of a charged pion. The structure of these two latter diagrams 
is different than the structure eq.(0.I) found for 7r° exchange, but it is easy to see that it is a 
combination of eqs.(01), (011) and after antisymmetrization with respect to the initial state leads 
again to the structures of eqs.(H). 

Should the KK system be in an intrinsic S-wave, L = 0, we would have eq.(|^II) and the cross 
section contains the factor p^, as correctly stated in ref. which largely affects the shape of 
the KK invariant mass distribution. 

As already mentioned, we will also consider tx^t] production. One interesting property of the 
TIT] system reflected in chiral dynamics is that it does not couple in P-waves to lowest order in 



the chiral counting [Q. It does not couple to vector mesons either ^J]. It can couple in higher 
orders but such effects are suppressed by more than one order of magnitude with respect to the 
dominant S-waves This means that there are no terms of the type of diagram 0c), d) or e) 



with the structure of eq.@, and in fact what one finds is that the matrix elements m^^^'* and u'^^^ 
of eq.(^ do not contain any vrvrr/ or r/r/vr terms. However, terms of the type of diagram |l]a), with 
the TTT] system in S-wave are allowed leading to a structure of the global amplitude of type II in 
eq.(||). The same comment applies to fig.|l|b) with respect to the K^K^ state. 

One still has to take into account that in the experiment only unpolarized observables are 
measured. This implies an equal probability for the initial state of being in any of the three 
possible total spin projections. Taking the vector Q to represent either — p^^o or p^, any of the 
structures shown in eq.(|P can be written as ^ = a^^^pi cr*-^-' It is straightforward 

to see that the matrix elements between states Ss) of well defined third component 5*3 and total 
spin 5 = 1 satisfy: 

(l/3|^|l«) = r^,^^|pi||Q|Yi„_^(Q) , (9) 

with Q the unit vector in the direction of Q, YLm{0, 4>) the usual spherical harmonics and the 
matrix rjajs is given by: 

/ +1 -1 \ 
r/= +1 -1 +1 . (10) 
V -1 +1 / 

In this matrix the rows correspond to a = +1, 0, —1, in order, and analogously for the columns. 
We have also taken pi parallel to the z-axis, that is, pi = (0, 0, \pi\). 

The structures and couplings of the amplitudes discussed here could be evaluated in explicit 
microscopic models. Given the large energies involved there could be many competing mechanisms 
some of which would require information, like partial decays of resonances, which is not available 
at present. Given these limitations our choice provides a reasonable and simple starting point 
from the phenomenological side. It also illustrates the use of Appendix A in order to take care of 
the final state interactions in any other model. 
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3 Final state interactions 



In this section we first discuss the FSI due to the meson-meson interactions and then we will also 
consider the FSI from the Kd channel. Afterwards, we sum up both contributions giving rise to 
our final renormalized amplitudes. 



3.1 Meson- meson final state interactions 



The K^K^ system in J = 1 will interact strongly and couple to the vr+r^ system. In the input 
of the lowest order chiral Lagrangian was used as the kernel (potential) of the Bethe-Salpeter 
equation which produced exact unitarization in coupled channels. The extension in [|T9| ,p4| to 
include effects from higher order Lagrangians, crossed channel dynamics, explicit exchange of 
genuine resonance states and scale independence did not modify the properties of the scalar sector 
found in |^ with only the lowest order Lagrangian and a cut-off of natural size, about 1 GeV, 
to regularize the loops. Since the divergences were only logarithmic the numerics are not changed 
when this cut-off is substituted by a regularization scale /z ~ 1 GeV and a subtraction constant 
that can be calculated in terms of the cut-off, see e.g. [|31],^,^ 0. Hence, given the simplicity 
and accuracy of ref. |2^, which is a limiting case of the more general formalism developed in 
(see e.g. [46, 4^), we will use this approach in our present problem. 

Diagrammatically it means that in addition to the tree-level diagrams of figjl| we will have the 
diagrams of fig.0 which contribute to the K~^K^ production in the first line and to tt'^t] production 
in the second one. By calling G the loop function of the mesons, the sums in fig.^] will dress the 
structure eq.(|^.II) containing the pd vector corresponding to the case when the two mesons are in 
S-wave in their CM reference system, L = 0. So we will have: 



V 



f^n\Pi\\Pd\'yimiPd 



\Pl\\Pd\'yim{Pd) {f^rj + fnrjG ~^ fKR^KK^KK^TTri) 



K+K'^ : f^j^lplWPdlYlmiPd) IPlWPdlYlmiPd) [JkR + f'KK^KKtKK^KK + f ^r^G ^^Kr,^ ^ k) , 

(11) 

where the subscript m corresponds to a- /3 in the notation of eq.(^, all the three- momenta refer to 
the CM frame of the pp system and fpq are the couplings of the two pseudoscalar meson systems 
with L = 0. The latter are defined such that the minus global sign in eq.(Bffl) is reabsorbed by 



them. Eqs.(ll) are more elegantly written in a 2x2 matrix form as: 



\pi\\Pd\YiUPd){l + tG) 



J IT T] 

■s 

KK 



(12) 



with tij the S-wave transition matrix KK, ttt] KK, vr?] in / = 1 with index "1" for tti] and 
index "2" for KK, and G a diagonal matrix G = diag(G'^^, G;^^). Now taking into account 
the Bethe-Salpeter equation and the on-shell factorization of the potential in the loop functions 



involved in the strong interactions proved in ||23|| : 

t = t2 + t2-G-t; t = [1 - t2 ■ G]~^- 12 



(13) 



"^While the explicit calculations have been done in dimensional regularization, this statement holds for any 
mass-independent regularization scheme. 
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where t2 contains the corresponding lowest order CHPT meson- meson scattering amphtudes |^ , 
we can write eq.(|12D as: 

J TT'r} 



\Pl\\Pd\yim{Pd) 



KK 



(14) 



where T> = [1 — 12 ■ G\. Both t and G are functions of the invariant mass of the meson- meson 
system, M|. The previous formahsm to take into account FSI was originally employed in the 
calculation of 77 — i>meson-meson in ref. ||2^ and later systematized in more general terms taking 
into account the analytic properties of the form factors in refs. ||2^,^. In these last references 
it can be seen that the above results for taking care of the FSI are exact when considering only 
the right hand or unitarity cut. We will take the functions f^^ and as constants since we 
are concerned in the energy region available to the ANKE collaboration which reduces to just 
about 45 MeV above the K^K^ threshold. Our approach of taking care only of the right hand 
cut corresponds to the expected dominance of the resonances ao(980) and A(1405) which are very 
close to the dK^K^ threshold. The chiral model of sec. p.2| gives rise to real couplings f^^ and 
and so we will take them in the following. Nevertheless, the results of Appendix A can be 



equally applied to real or complex coupling functions /^^ and 



KK- 



3.2 Anti-kaon-deuteron final state interactions 

Now we consider the FSI from the Kd system. The interaction of the with the protons 
and neutrons is rather weak and we will neglect it. However, this is not the case for the 
K^n interactions which are very strong close to threshold due to the A(1405) resonance below 
the K^n threshold P5|-p2|. On the other hand what we need here is the interaction with 
the deuteron that is quite strong close to the threshold due to extra reinforcement of the multiple 
scattering of the K in the deuteron as proved in multiple evaluations of this quantity using Faddeev 
equations P9|-|53[] . A reanalysis of this quantity to the light of the new KN amplitudes generated 
in the chiral dynamical approach of was done in within the fixed scatterer approximation 
for the deuteron, which proves rather accurate comparing the results with those of the non-static 
calculation of 



A sizeable Kd scattering length of about (— 1.6 + zl.9) fm is obtained in [Q. In 
order to take into account this extra interaction we first extrapolate the results of the Kd scattering 
amplitude at threshold of |Q to the small finite K energies of the ANKE experiment [||]. For this 
purpose we rely upon the results for the KN scattering matrix found in [^(fig.9) which show 
a drastic reduction of the real part of tx-p + tK-n at ^/s ~ 1450 MeV. Taking into account this 
fact plus the approximate good results of the impulse approximation for the imaginary part of the 
Kd scattering length and the fast decline of the real part of the amplitudes from threshold on, 
suggest a quadratic interpolation between the results at threshold and the impulse approximation 
at a/s — 1450 MeV and beyond. Hence, for the general and illustrative purposes of the present 
chiral model to the primary production amplitudes of the KK and ttt] channels, the following 
parameterization is expected to provide a sufficiently accurate description of the Kd scattering 
amplitude: ^ ^ ^ 

Re tKd{MB) = a{MB - MboY + KMb - Mb,) + c , (15) 
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with 



a = 4.32 ■ 10-^ MeV"^ , b = -1.55 ■ 10"^ MeV'^ , 
c = 0.13 MeV~^ , Mbo = 1432 MeV , 



(16) 



where = {p^o — is the invariant mass of the K'^ and the neutron in the deuteron. The 

imaginary parts are well approximated by the impulse approximation and are almost constants in 
the whole interval. We take finally: 

Im txd(MB) = b'{MB - Mb,) + c' (17) 

with 

6' = 1.1 ■ 10-3MeV-2 , c' = -1.5- 10-1 MeV"^ . (18) 
We apply the formulae ([T|) and (|1^ for Mb < 1.45 GeV, and 

tUMB) = t^d(1.45 GeV) , (19) 

for Mb > 1.45 GeV. We have also checked that substituting this limiting value by a larger one 
only modifies mildly the resulting distributions. 

The implementation of the FSI of Kd requires to rewrite the amplitudes of eq. (^ , for Q = 'Pk+ — 
Pj^o and Q = Pain the rest frame of the Kd. This is easily done by taking into account momentum 
conservation and the fact that pk+ —Pk'^ is a Galilean invariant. Given the small velocities involved 
in the dK^K^ system we find it appropriate to just apply Galilean transformations. Let us first 
consider the intrinsic P-wave meson-meson contribution, eq.(H.I). In the following we denote by 

= Vk+ —PRo, and hence this contribution involves YimipG-)- Since we are considering Galilean 
invariance p^ = p'_^, we then can rewrite |p^|Yim(p^) as: 



|p^|Yi„(p^) = |p^|Yi„(p'^) = |p/^+|Yi„(p'^+) - |p'^o|Yi„(p'^o) , (20) 

where the primes stand for variables in the K^d rest frame. The term YimiP'j^o) involves the 
P-wave contribution of the K^ and hence we neglect its modification in the present case of low K^ 
energies. The term Yim{p'j^+) does not depend on p'^o- This implies that the deuteron and K'^ are 
in a relative S-wave and can suffer the strong Kd interaction. The Kd FSI are diagrammatically 
represented in fig.^ and the corresponding terms are renormalized by changing them by: 

l + GdtRd, (21) 

where G^ is the meson-deuteron loop function for the KN interaction. However, although the 
deuteron effects appear only as recoil corrections with respect to the nucleon meson-baryon Gn 
loop function these effects can be around mx/Md — 25% due to the large mass of the kaon, cf. 
fig.^D). The function is given in |^ using dispersion relations and a subtraction constant 
'^KN = —1.82 as needed in the approach of to reproduce the low energy KN results of pO] 



using a cut-off to regularize the loops. Because of the already mentioned recoil effects and the 
large momentum transfer (through the shaded areas in fig.0), we identify Gd = Gn but allow for 
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a variation of ~ 30% in a^N- Therefore, we present results for a^^ = 
said this, we find that eq.(^) is renormahzed by the FSI as: 

Now, taking into account the following equalities: 



-1.84 and —1.3. Having 
(22) 



Pk+ 

Pro 



Pk+' 



Md + rriK 



PRO + Pk+ 



ntK 



Md + niK 



(23) 



it is then straightforward to rewrite eq.(p2D as: 



fKR\Pl\ 



\PK+\YlmiPK+) 2 + -GdtRd + \Pd\Y ImVPd) 

Md + niK 



(24) 



where we have multiplied the previous structure by the K^K^ P-wave coupling times \pi\-, 
since both factors appear in the original production process. In addition, Md = 1875.61 MeV 
is the deuteron mass and rriK = 495.7 MeV = {mK+ + mj^o)/2, with mK+ = 493.677 MeV and 
rrij^o = 497.672 MeV the and masses, respectively. 

We can follow similar steps to take into account the FSI of the K^d to the CM meson-meson 
S-wave contribution, L = 0, eq.(||.II). Considering the identity: 



Pd = -Pk+ 



Md 



it then follows: 



\Pd\Ylm{Pd) = - \PK+\'y{p'K 



Md + 2mK 
Md 



Pro 



\p'ro\Y{p'ro) . 



(25) 



(26) 



Md + rriK 

As discussed above the Y(p'^o) term involves pure P-wave and it is not renormalized by the strong 
K^d interaction. Its full contribution is then accounted for by eq. (p^ . Hence only the term 
proportional to Y(p'^+) = Y(px+), eq.(p3D, is renormalized due to the S-wave K^d interaction as: 



- fKRWPlWPK^mpK^] 



Md 



-Gdt 



d l^Kd 



Md + rriK 

multiplied by the corresponding coupling constant already introduced in eq. (p!l|) . 



(27) 



3.3 Renormalized amplitudes 



Once we have taken into account the important FSI due to the resonant meson-meson and K^d 
interactions, eqs.(p^), (^) and (^), the renormalized dn^rj, F^^+n, and dK~^K^, F^+^o, production 
amplitudes, corresponding to the transition between total spin third components a ^ f3, read: 



4v/^' 



\pi\\Pd\Yic.-p{Pd) [D-\M^)]ufl + [V-\M^)]ul 



eS 

' KK 
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40F 



+ Val3\Pl\\PK+\'yia^l3{PK 

Md + 2mK 



Md + rriK 



Md + niK 
Gd{Ml)tUMB) 



f^^Gd{Ml)tfid{MB) 



(28) 



The double invariant mass distributions are obtained straightforwardly from the previous equa- 
tions after summing over the final state polarizations and averaging over the initial ones. In this 
way one has: 



= lQnC^-^,e{l-\cose^,\)MjMs\pd\' \[V-\M^)U^^ + [V-\M])],,l 



S |2 

kk\ ' 



16nC^-^9il-\cos9K+\)MiMB\ \pd\' \P'\M^)]2if^^ + [V-\M^)]22f^j^ + 

Md + 2mK 



dMjdMB s3/2 



+ \P 



P |2 

kk\ 



K+\ 



Md + rriK 



2 + 



Md + rriK 

+ 2|p,||p-W+|cos^K+ Re ([V-\Mj)Uf^^ + [V~\Mj)]22f^^ + fKK 



Gdm)tKd{MB) 



-Md 



Md + rriK 



fKKGd{Ml)t^d{MB) + f, 



p _ 

KK 



2 + §^G,(M|)*,,(M«) 

Md + niK 



(29) 



with s = {pi + P2)'^ for the two protons in the initial state {^/s = 2912.88 MeV for the ANKE 
kinematics considered here), Mb is the corresponding rjd or K^d invariant mass and C is a nor- 
malization constant. The cosine of the angle between Pt,+ {pk+) and pd^ cos 6',,+ (cos 6*;^+), can be 
written in terms of M/, Mb- By integrating with respect to M/ or Mb in eq.(p9D we can obtain 
the invariant mass distributions with respect to Mb and M/, respectively. 



4 Results and discussion 

Apart from the absolute normalization of the amplitudes, our chiral model for the primary pro- 
duction depends on two free parameters, 9 and 0, such that 

ffik = cos9 , f^^ = sin 6* cos , /^^ = sin 61 sin . (30) 

First, in order to show the relevance of the FSI, we take 9 = (p = 0, implying f^^ = and 
J^KK ~ ^- fig-0 we display several curves corresponding to dax+^o/dMj neglecting either the 
two considered FSI and including either one or two of them. In the rest of this section we take 
ci'KN = —1.84 unless the contrary is explicitly stated. The distribution in the absence of any FSI 
(dotted line) peaks around M/ = 1003 MeV. If the K~^K^ interaction is switched on (dashed 
line) the strength is shifted considerably towards low invariant mass and the peak moves to about 
Mj = 997 MeV. This is an obvious consequence of the presence of the ao(980) resonance around 
980 MeV and the K+K^ distribution feels the tail of that resonance which increases the strength 
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the closer one is to the resonance position, and hence to smaller values of the K^K^ invariant 
mass. If one switches on only the K^d FSI (dashed-dotted line) the distribution is rather broad and 
there is an accumulation of strength to higher values of the M/. Finally, when all the interactions 
are considered (thick solid line) the peak of the distribution moves back to lower masses around 1 
GeV where the pure phase space peaks as well. The strength is furthermore increased by about a 
factor five due to the combined effects of both FSI. The effect of the K^d interaction moving the 
peak towards the center of the distribution reflects the fact that for these values of Mj the 
and the deuteron are at rest where Gdixd has its maximum. Indeed, in the extremes of the M/ 
distribution either the kaons go together and the deuteron goes opposite to them, or the deuteron 
is produced at rest and the two kaons go back to back. In both cases the K^d invariant mass 
is relatively far from the K^d threshold situation. In addition, we also include the full result for 
(^KN = —1.3 presented by the thin solid line. As it is clear this variation in the value of aj^^ 
mostly decreases the width of the distribution while the peak position is just slightly decreased by 
less than 2 MeV. Had we further reduced the value of a^jv; instead of increasing it, the changes 
would have opposite sense. 

Similar changes in strength and shape can be seen in flg.^for the dax+Ko/dMs invariant mass 
distribution where the notation for the curves is the same as for flg.^ It is remarkable that 
the strong K^d interaction at threshold pushes the mass distribution towards lower Kd invariant 
masses as a reflection of the presence of the A (1405) resonance in the K^n system, much as 
in the case of the K~^K^ distribution where the presence of the ao(980) resonance pushed the 
distribution towards low K^K^ invariant masses. Here the effects of decreasing the modulus 
of opposite to those in the K^K^ mass distribution pushing the distribution to higher 

invariant masses. 

In flg.|^ we show as a function of 6 and the peak position of the daj^+j^o/dMj distribution 
and its width. The width is deflned to be the difference between the values of Mj for which the 
number of events is half of the maximum value. The same is shown in flg.0 for a different value 
of the subtraction constant aj^j^, i.e. aj^j^ = —1.3. The advantage of these figures is that one can 
approximately describe the shape of the distribution as a function of 6 and (j). 

Next, we investigate the role of the 6 and parameters in the total production of tt^i] and 
K~^K^. The TT"*"// production is mostly done around the ao(980) resonance region. In fig.^] we 
show the ratio between the integrated it^i] production cross section between Mj = 950 MeV 
and the end of its phase space and the K^K^ production cross section in all its available phase 
space. We can see in the figure that for most of the values of and the vr"*"?] production rate 
is substantially larger than that of K^K^. It is interesting to point out that even in the case 
when there is no primary vr+r] production so that f^^ = (6' = and any or = 7r/2, 37r/2 
and any 6'), the final state interactions starting from primary K^K^ production can lead to a Tt^r] 
cross section an order of magnitude bigger than that of the K^K^ . One also finds interesting 
interference effects for some values of 9 and that can reinforce the n^r] production as compared 
to the K^K^ as well as other situations when the K^K^ is produced more copiously than the tt^t] 
channel. The former occurs for values of around vr and 2tt and of 6 around 7r/2 so that ~ 
and f^j^ ~ with o"(7r+?7) > ?)Qa{K^K^), while the latter tends to happen for a wide range 
of parameters leading to K^K^ invariant mass distributions peaked at values of Mj higher than 
1010 MeV. This spectacular dependence of the ratio of ti^t] to K^K^ production on the values of 
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primary production weights 6 and should obviously serve as a stimulation for the experimental 
measurement of the vr+r] production cross section. 

It is interesting to see also the shape of the vr+r^ invariant mass distribution, which is independent 

of 

(^KN- We show in fig.^by the dashed-line that the normalized daT^+jj/dMj event distribution for 
f^^ = and f^j^ = arbitrary (|/^^| < 1, see eq.(^0|)) has no clear signal of the ao(980) resonance 
around the values of Mj = 980 MeV. This seems somewhat surprising since the coupled channel 
approach definitely generates the resonance and we have already observed the effects of its tail 
in the K~^K^ invariant mass distribution. The lack of resonance structure is due to the P-wave 
character of the reaction and the appearance of the Ip^p factor in This factor grows as 

the invariant mass decreases and distorts the ao(980) shape. In fact it is interesting to observe 
that if we divide (ia^+^/rfM/ by |pczp, which is also shown in the figure by the thin solid line, the 
resonance shape appears with a width of around 40 MeV and with the peak at 990 MeV. Noticing 
this fact is also important from the experimental point of view in order to extract properties of 
the ao(980) resonance in this reaction. 

We can see that the distributions are rather dependent on the values of the 6 and 0. This fact 
could be used to extract the optimal parameters from the data on K^K^ distributions, assuming 
that a good fit is possible. Should this be the case, the theory would then predict absolute rates 
and mass distribution for the tt+t] production or other experimental yields, which would be a real 
prediction of the approach in spite of having started from two unknown parameters. 

5 Conclusions 

In this paper we have performed a phenomenological study of the pp —>■ dK^K'^ and pp dn^r] 
reactions close to threshold presently studied by the ANKE collaboration at COSY. We have 
emphasized the relevance of the final state interactions which is quite important in the present 
case due to the proximity of the K~^K^ system to the ao(980) resonance and the K^n system 
to the A (1405) resonance. We found that the consideration of these interactions has important 
consequences both in the shape and strength of the invariant mass distributions. We also studied 
the interaction of the two final states K~^K^ and n^i] by means of a coupled channel chiral 
unitary approach which generates both the ao(980) and A(1405) resonances. Given the freedom 
in the primary production amplitudes we parameterize them in terms of three types of structures 
involving i = 1, L = and i = 0, L = 1 for the K~^K^ channel and the only allowed i = 1, L = 1 
for iT~^rj production. This left us with two independent parameters (up to a global normalization 
of one cross section) together with a subtraction constant with an expected uncertainty of about 
25%. The sensitivity of the shapes of the K^K^ and K^d invariant mass distributions to those 
parameters was investigated anticipating that the measurements of these quantities could serve to 
fix them. This would allow us to make absolute predictions for tt^t] production due to the dynamics 
of coupled channels generated in the chiral unitary scheme followed here. Other quantities which 
might be measured could also be predicted in that case. Furthermore, we observed that the tx^t] 
production was dominated by the ao(980) resonance and that a clear signal for the relevance of the 
K^d FSI would be the observation of a peak towards low K^d invariant masses in the da^+K^ / dMs 
differential cross section. On other hand, we have also pointed out that the ao(980) would not be 
clearly visible in the data for da.„+n/dMj because of the \pd\'^ factor due to the P-wave character 
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of the reaction which distorts the shape of the resonance. Yet we found that the shape of the 
resonance was regained by dividing da.„+r^/ dMj by 

Finally, we have also provided general expressions to take into account the FSI derived in this 
paper for any other more specific model of the primary production mechanism. 

The study done here clearly shows how the measurements performed or planned with those 
reactions provide basic information on the strong interaction underlying the meson-meson and 
meson-baryon dynamics and should produce complementary and valuable information to the one 
obtained from other processes. 
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A General structure of the process pp — > dPQ 

Let us denote by i the relative orbital angular momentum of the deuteron and the CM motion of 
the two pseudoscalar system PQ and by L the orbital angular momentum of the latter in their 
own CM frame. As discussed in sec. ^ close to the dK^K^ threshold the leading contribution 
stems from £ = 1, L = and £ = 0, L = 1. We denote by ^Q the orbital angular momentum of the 
two protons in the initial state and, as also discussed above, the only possibilities are £o = 1? 3. 
Finally, the symbol S'^ refers to the total spin of the deuteron, S''' = 1, with its third component 
indicated by S^. Analogously 5* refers to the total spin of the pp system, which is also fixed to be 
1, and 5*3 indicates its third component. 

Keeping only the components relevant for our reaction, we can consider the following angular 
momentum decomposition of the final dPQ state as: 

\d{Sipa)P{k)Q{k2)) cc Y,C{S;^m\llJ)(Y,UPdr\J,S3 + m-i=l,L = 0) 

m,J 

+ Yi^(fc)*|J,53 + m;£ = 0,L= 1)) +... (A.l) 

where k is the PQ-CM three-momentum of the pseudoscalar P, the symbol C{mi m2\ji 32 J) is 
the Clebsch-Gordan coefficient for the composition of two angular momenta ji and j2 to give the 
total one J and the ellipses simply denote other terms of no interest here. It is worth noting once 
again that for the n^rj system only the ^ = 1, L = component is relevant due to the absence of 
resonant interactions of this system with the deuteron. Performing an analogous decomposition 
for the initial pp state we can write for the transition matrix element: 

{d{SiJd)P{k)Q{k2)\T\p{p^)p{p2),S:,) = J2 CiS^ S3 - 53^11 1 J) 0(83 0|1 io J) Yio(pi)* 

J,£o 
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X Y 



(A.2) 



where pi = (0, 0, \pi\) and because of the small velocities involved at around the dK^K^ threshold 
we can also write Yim(p^) instead of Yim(/c) with defined in sec. ^. Note as well that we take the 
invariant matrix elements T^-^^^ = 0, due to the absence of any resonant S-wave interaction between 
the 7r~^rj system and the deuteron. In order to take care of the final state interactions due to both 
the meson-meson and Kd S-wave interactions one can proceed in a completely analogous way to 
that of sec. 0. First we define the related quantitides T^qi^ = \Pd\^ioe^ and T^^^^ = \p^\Aq^^ 
and second the A'^^^ can be taken, if desired, as constants. In this way one has: 

= \pd\Y^s,-si{Pd)Y,Y,,,{p,r C{Si S, - Si\l 1 J) C(53 0|1 £o J){[V-\M^,)]nA 



Jtttj 
10£o 



+ [V-\M',)UAi,l^) , 

= \mis,-si{Pd)Y.^upiy c{si s, - si\i i j) c{s, o\i 4 j)([v-\m])],,a 



j,eo 



iJ-rrrj 

hoeo 



+ [V-\M])h,A(,l^ + Al^^^ 



xCiSs 0|1 io J) 
Md + 2mK 



OUo J 



Md + rriK 



+ \PK^\yis,-siiPK+)J2^ioo{Pir c{st s, - s',\i 1 J) 

GdiMl)tf,d{Ms)Kl^ (A.3) 



+ 2 + 



Md + rriK 



GdiM],)tKd{MB)]A^ 



\JKK 



Once a model for the primary production mechanism of the PQ systems is developed, the 
functions T^j^^ can be determined and from them the FSI state interactions can be taken into 
account by eg. ( |A.^ ). We can do this exercise for our previous model by comparing eq.( |A.2| ) of 
the present appendix with eq.(^, times the appropriate couplings constant fpq or /^^. Taking 
£0 = 1, since our model only involves one power of pi, we then have: 



IKK _ Sy/^ fP 



rp2KK 

101 







rp2 KK 
-^010 



(A.4) 



Twi = -4^V^fnv\Pi\\Pd 
. 



^ 101 



Tiri = '-§f!,mpd\. 
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Figure 1: Chiral model for the primary production used to extract the structures given in eq 
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Figure 2: Diagrams relevant to take into account the meson-meson FSI. The ao(980) resonance is 
dynamically generated through the iteration of the meson-meson bubbles. This iteration is indicated in 
the figure by the ellipses. 
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Figure 3: Diagrams to take into account the FSI. 
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Figure 4: da{K~^K^)/Mj for the whole range of available M/ in the reaction pp —>■ dK^K^ with 
^ = 2912.88 MeV. The thick (thin) sohd hne is the full result with a^j^ = -1.84 (-1.34). The dashed 
line corresponds to including only meson-meson FSI, the dashed-dotted one includes only K^d FSI and 
the dotted line includes no FSI with a p^ factor for the modulus squared of the amplitude. 
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Figure 5: da{K^ K^)/Mb for the whole range of available Mb in the reaction pp dK~^K^ with 
^/s = 2912.88 MeV. The thick (thin) soHd Hne is the full result with a^N = -1-84 (-1.34). The dashed 
line corresponds to including only meson-meson FSI, the dashed-dotted one includes only K^d FSI and 
the dotted line includes no FSI with a factor for the modulus squared of the amplitude. 
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Peak Position in GeV 





Figure 6: The upper panel shows the peak position of da{K'^ K^) / dMj and the lower one the width of 
the same distribution for a^j^ = —1-84 as functions of 9 and (f) (in radians). 
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Figure 7: The upper panel shows the peak position of da{K'^ K^) / dMj and the lower one the width of 
the same distribution for O jy — 1.3 clS functions of and (j) (in radians). 
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Figure 8: a{TT'^ri)/a{K'^K^) as function of 6 and (f) (in radians). 
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Figure 9: da{Tr~^7])/dM[. Solid line, full result, dotted line does not include FSI, with a factor from 
the modulus squared of the amplitude. The thin solid line corresponds to the full result but divided by 
times 250^ MeV^ (to normalize the curve to the full result at the KK threshold). 
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